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Abstract 

The letter deals with the spin-freezing process observed by means of NMR- 
NQR relaxation or by muon spin rotation in underdoped cuprate supercon- 
ductors. This phenomenon, sometimes referred as coexistence of antiferro- 
magnetic and superconducting order parameters, is generally thought to re- 
sult from randomly distributed magnetic moments related to charge inhomo- 
geneities (possibly stripes) which exhibit slowing down of their fluctuations 
on cooling below Tc . Instead, we describe the experimental findings as due to 
fluctuating, vortex-antivortex, orbital currents state coexisting with d-wave 
superconducting state. A direct explanation of the experimental results, in 
underdoped Yi_^Caa;Ba2Cu306.i and La2-2:Sra;Cu04, is thus given in terms 

of freezing of orbital current fluctuations. 
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NMR-NQR and /iSR experiments show that in underdoped cuprates, on coohng from 
about Tc, divergent behaviour of the relaxation rates occur. This phenomenon, generally in- 
terpreted in terms of "glassy spin- freezing" or of "coexistence of superconductivity and mag- 
netic ordering" , is believed to be related to magnetic moments resulting from charge inhomo- 
geneties, possibly stripes. For ^^^La NMR-NQR observations in La2_a;Sr2.Cu04, see Julien et 
a/.0 and References therein; a review of early data is given in Ref. ^; the muon longitudinal 
relaxation rate and spin precessional frequencies in La2-:rSr^Cu04 and in Yi_xCaxBa2Cu306 
have been measured by Niedermayer et a/.i. More recently the low-temperature spin-freezing 
process in underdoped Yi_2;Ca2;Ba2Cu306.i ( Tc ~ 35 K) has been studied by "^^Y NMR 
relaxation in a field Hq = 9.4 T(Ref. |). The occurrence of a fast-relaxing component in 
the ^^Y spin-lattice relaxation below T^ confirmed the /zSR findings! interpreted in terms 
of coexistence of superconductivity and spin-freezing process. A new significant aspect was 
the detection of the stretched character of the relaxation timeS, while the /iSR data! had 
been discussed in terms of a single exponential recovery. 

Very recently Mook et a/.i by means of inelastic neutron scattering measurements in 
underdoped superconducting YBCO detected longitudinal with respest to c-axis magnetic 
moments of unknown origin , with fluctuation frequencies well below the energy resolution 
of 1 meV. Similar observation, but with transverse magnetic moments , was reported by 
Sidis et a/.i. NMR-NQR and /xSR spin-lattice relaxation dataBH^ have been interpreted in 
terms of fluctuations of a local h{t) originating from magnetic moments associated to hole 
localizationi'ior from staggered moments within locally ordered mesoscopic domains!, with 
non-zero effective spinlli or from stripesS. The microscopic origin of the fluctuating field h{t) 
is still an open issue. Here we address this problem within a substantially different picture 
, namely the one of the extended charge density waves with imaginary order parameter 
{id-CDW) . Important features of this new scenario are described here and it is pointed out 
that the fluctuating field originates from vortex- antivortex orbital currents. 

We start from the so- called t — J model by including the inter-site Coulomb interaction: 
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with 5i being the number of extra holes per unit cell of bilayer, ipf'^''^ , [ipj''^^) creation 
(annihilation) operator constructed on the basis of the singlet combination of copper {d) 
and oxygen states (p) ( see Ref.§ for details). It should be noted that the temperature 
dependence of Cu(2) Knight shift in YBa2Cu408, the evolution of the Fermi surface in 
Bi2Sr2CaCu208+j; and the k-dependence of the pseudogap in the normal state have been 
successfully explainedlli in the framework of that model, including the d-wave orbital sym- 
metry of the superconducting state (SC)0. The d-wave SC can coexist with extended charge 
density waves (for short id- CDW)li@. Chakravarty et a/.li3 have stressed that the corre- 
sponding ordered phase is a staggered pattern of orbital currents. The region in which d-wave 
SC coexists with id-CDW strongly depends on dopingS. For Sopt the SC phase forces out 
zd-CDW state (see Fig. 1 in Ref. H), but when T*/T^ <2.5 then SC and id-CBW state can 
coexist even at low temperatures. Non-monotonic temperature dependence of energy gap 
parameter, due to the competition of SC and id-CDW, was predicted in Ref.lli and confirmed 
by Ekino et a/.lli by means of break-junction tunneling spectroscopy in Bi2Sr2CaCu208+y. 
Starting from another form of Hamiltonian Murakamilll has predicted (for G < 0) the coex- 
istence id-CDW, antiferromagnetic (AF) and r^-triplet states near half filling. Coexistence 
of d-CDW and spin density waves (SDW) was also derived by Bonis et a/.0. 

The main issue related to the scenario recalled above is whether it is possible to have a 
coexistence of d-SC and AF or SDW in the frame of the singlet - correlated band Hamiltonian 
(|^). Using the Roth - type decoupling scheme^, in a mean field approximation, from Eq.(l) 
one has 

^MF ^ ^ M;+M,v^, (2) 

k 

where the matrix Mu is 
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Vk+Q-Q ^k+Q — 

\-iU^k^Q,QY -(Aifc.g)* -W-k-Q,QY -S^k-Q + l^) 

Let us focus now at the order parameters: 
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where G^fc are the Fourier transforms of transfer, superexchange and Coulomb couphng 

parameters, respectively, eg and Sq are the Fourier amplitudes of the conventional charge 
density waves and spin density waves. For the commensurate instability wave vector Q 
= (7r, tt) as it was observed in neutron scatteringii one derives the following relations 

^I,Q = ^Ik + ^2k + ^-^I 

vl+Q-Q = ~iSlf.)* + {Sli, + iDl)* 
iV-k-q)* = ^ik + ^2k ~ '^P^k 

W-k-Q,Qr = -{si,r-{si,-iDir 



(6) 



Let us consider now the quantities 
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where the indexes 1 and 2 refer to the first and the second nearest neighbour hopping param- 
eters. Slj. changes the sign under the transformation k — > k + Q and as one can see from 
Eq. this yields a damping factor, since the energy dispersion becomes complex. The 
value 5*2/; is invariant under such transformation and can generate s-wave CDW and SDW. 
Therefore, for open Fermi surface one finds that at the hopping parameter ti larger enough 
the expectation values < Cg > and < Sq > must be zero, and conventional SDW/CDW do 
not exist. This conclusion is in agreement with the one derived in Ref.il on the basis of 
neutron scattering data. However, we stress that the above conclusion is valid only for the 
energy dispersion of the form (cos + cos ky) , but not for the case C/^ — (cos k^ ~\~ cos ky^ 
when pockets are formed around the points (±7r/2, ±7r/2) in the Brillouin zone. In our 
picture this happens when antiferromagnetic correlation (parameter J) hampers the nearest 
neighbour hopping i. e. ti [Ppd + %~) ~ 0- It is noted that the order parameter Uk resem- 
bles the so-called //-singlet pairing introduced by Yan^. Formally the f/^-order parameter 
construction in Eq.(^ also resembles the Larkin-Ovchinnikov, Fulde-Ferrell parametei^iii 
with the important difference that the instability wave vector q is completely different. The 
leading term of Uk 

ul = -j^T. + ty) (^I;*^i'f^Q> (9) 

^pd k' 

is determined by the second- and third- neighbours hopping integrals. Numerical solution 
of Eq. (^ in correspondence to the parameters (^1=80 meV, ^2 = and ^3= 12 meV) usecS 
for the Fermi surface near optimal doping, yields a critical temperature of the onset U -phase 
below 1 K. One should remark that the onset of this inhomogeneous superconducting state is 
very sensitive to magnetic impurities and to the details of the Abrikosov's vortex lattice. In 
addition we note that the ^/-ordering is particle-particle pairing at opposite spin orientations 
and must be sensitive to the external magnetic field. Since there is no evidence of such an 
effect of the magnetic field (see the data later on) we conclude that the peak in relaxation 
rate is not related to phase transition involving this U order parameter. In the matrix given 
by Eq. (^) there are two imaginary components iD"^ and iD^ It is useful to introduce their 
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combinations 



the first corresponding to charge current and the second one to spin currentil'ii or, in 
other terminology, to the spin-nematic stateil'ii. It is clear from Eqs. (|TO|)-(|Tl|) that 
T^^(ic/- CDW) » T'='^(irf-SDW), and when 2G ^ J Eq. (|ll]) has no solution at all. 
Concluding this theoretical analysis, we see that the real situation is strongly dependent 
on the Fermi surface. Because there are no indication about pockets-like Fermi-surface in 
superconducting underdoped cuprates, we do expect coexistence of the d-SC, zd-CDW and 
id-SDW or fZ-state at low temperatures. Let us see now what the experimental findings tell 
us . 

On general physical grounds NMR-NQR or muon relaxation rates can be written 

l/T,^jl{hlff)j{u;m,Te) (12) 

where (^e//) is a mean square amplitude of the transverse effective field at the nuclear or 
muon site and (cj^, Tg) is the spectral density of the fluctuations at the measuring frequency 
ujrn and an effective correlation time Te is assumed. The stretched exponential recovery, of 
the form exp [— t/Tf]^''^ observed experimentallyS'l is naturally explained by charge density 
waves scenario, as it was pointed out by PhilipJ^. Usually this process is described as a 
superposition of Debye relaxation or in an ideal limit as a Laplace transform 

1/2 



exp 



t 

'Tf 

1/2 



/ p(l/Ti)expH/Ti]d(l/Ti) (13) 
Jo 



where p(l/Ti) = ^ (^) exp 



Ti 
4Tf 



is the distribution function. It is conceivable to 
assume that the distribution of (1/Ti) is related to a distribution of correlation time and 
in turn to a distribution of the energy barrier E pinning the fluctuating - sliding current 
motions. Then we write 



{T,) = ToeM{E) IkBT) (14) 

and 

J{u:^.r;)= ^^l^] (15) 

1 + '^m (re) 

When the progressive slowing-down , on coohng, of the orbital current excitations causes 
< Tg >~ cu"^, one has the maximum in the relaxation rates (frequently used to define Tg, 
commonly called " spin -glass" freezing temperature). From this maximum one can extract 
the effective magnetic field induced by the currents at the nuclear or muon site: 

1/2 



In 



(Tf 



(16) 



One could remark that on the basis of the above equation a slight underestimate of y (^e// 
is obtained, since the distribution of correlation times tends to level the maximum of the 
spectral function. The distribution also accounts from the deviation of the experimental 
data from the simple law ( Eq. ([T^)) in the low temperature region (see Fig. 1). The values 



of Y\^e/// derived from the experimental data on the basis of Eq. (|T^) are collected in 
Table 1. 

In Fig. one sees that Eq. (|12]), with Eqs. (p!^ -(p!5[) rather well fits the observed tem- 
perature dependence of relaxation time in spite of the roughness of the assumptions. The 
energy barrier < ii^ > is roughly the same, for different kind of experiments and systems (see 
also Table 1). It is naturally related to periodic lattice potential. It should be stressed that 
the energy barrier is almost independent on the applied magnetic field. This fact would be 
hard to explain in the spin-freezing scenario. On the contrary, in our picture the correlation 
functions entering Eq. (|10|) are of particle - hole type and they are proportional to difference 
of Fermi functions with the same spin orientationEl. Therefore, Zeeman energy comes out 
from these functions and they become independent of the external magnetic field. As a re- 
sult the pseudo- gap temperature T* and the pinning process caused by the lattice periodic 
potential are insensitive to the magnetic field. 



7 



As one can see from Table that the order of magnitude y \^lff/ of ^^e magnetic field 
agrees with theoretical estimation made on the basis of similar bond currents patterns (flux- 



The z— axis of the electric field gradient at ^^^La site is almost parallel to c— axis ( see 



for Ref. @). Therefore we can interpret ^(^hlj^j =255Gs as a field transeverse with respect 
to c-axis. The effective field at Y site (9Gs) is relatively small. This means that transeverse 
components within Cu — bilayer are almost compensated. Two magnetic field components 
at nuclear spin occur: a direct component - h{l) as dicussed above and so-called hyperfine 
and supertransferred fields from adjacent Cu(2) states - h{2). The ici-CDW state contributes 
to both. The modulation of the hyperfine coupling yields the effective magnetic field which 
is proportional to charge-spin correlation function 



where F[q) is the form factor involved in the Mila-Rice Hamiltonian (see Ref.H). This factor 
is different for La and Y nuclei sites and for electronic shell of ions Yly^^ and Er'^^ as well. 

In summary , we have described the so- called spin-freezing process in superconducting 
state of underdoped cuprates with non-zero first neighbours hopping in a frame of a regime 
of coexistence of d— SC and id— CDW states . 

We have shown: 

i) that this scenario is naturally derived in a frame of singlet - correlated band model 
(or on the basis of the t-J model including inter -site Coulomb repulsion) and it is consis- 
tent with recent neutron scattering data and tunneling spectroscopy measurements, ii) the 
progressive slowing down of the sliding current motions, due to pinning barrier related to 
periodic lattice potential explains the temperature behavior of the relaxation rates detected 
in NMR-NQR, /iSR and EPR experiments; in particular it directly justifies the stretched 
exponential character of the recovery curves; the fluctuating frequency lying in 10-1000 MHz 
range explain why the resolution limit (~ meV) makes them invisible in neutron scattering 
measurements iii) the fact that practically no difference is observed in the relaxation rates 



phase) Ref.ii. 




h{2) oc F{q - q)eqSq' 



(17) 
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upon increasing the magnetic field from zero (NQR) up to 23 Tesla and when the field is 
applied parallel or perpendicular to the c- axis ( that is not explained in a spin-freezing sce- 
nario ) is explained by the insensitivity of the current correlation functions to the external 
magnetic field . 

Thus we believe that our picture brings a new and suggestive perspective, which is fully 
supported by the experimental findings and reveals new insights on a much debated issue. 

Useful discussions with R. Zeyher, P. Carretta, M. Corti , A. Dooglav ,1. Eremin and 
F. Tedoldi are acknowledged. M. V. Eremin is grateful to Landau-Network-Centro - Volta 
foundation for the fellowship and possibility to stay at University of Pavia.His work in Russia 
was supported by Swiss National Science Foundation (Grant N 7SUPJ062258) and Russan 
Program " Superconductivity" (Project N 98014-1). The work is a part of the project of 
advanced recearch ( PRA - SPIS) of INFM. 
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FIGURES 
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1/T (K"') 

FIG. 1. Temperature behavior of the relaxation rates in La2-xSr2;Cu04 and in 

Yi.85Cao.i5Ba2Cu306.i, normahzed to the values at T>>Tg ?a Tmax, as a function of T^^. The 
data are taken from Ref.S'i: circles are for Lai.94Sro.o6Cu04 in zero field and resonance frequency 
18.6 MHz (NQR), stars - Lai.9Sro.iCu04 in H=23.2 Tesla, parallel to the c-axis, uptriangles - 
Lai.9Sro.iCu04 with H=9.4 Tesla, perpendicular to the c-axis, downtriangles - Lai.9Sro.iCu04 at 
H=9.4 Tesla, parallel to the c-axis, squares - Yi.85Cao.i5Ba2Cu306.i in H=9.4 Tesla. 
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TABLES 



Table 1 Maximum values of the relaxation rates, estimated values for the effective field and 
for the average energy barrier for pinning of the sliding motions of the orbital currents. 



sample 


Method (Ref.) 


I/O (MHz) 


Ho(Tesla) 


l/Ti(s-i) 


< E >(K) 


y< hi >(Gs) 


Lai.9Sro.iCu04 


i39La NMR0 


139.528 


23.2 


^700 


19 


210 


Lai.94Sro.o6Cu04 


i39La NQR§ 


18.6 





«8000 


30 


255 


Yo.85Cao.i5Ba2Cu306.i 


nmrI 


19.61 


9.4 


ssl 


26 


9 


Yi„^.Cax.Ba2Cu306.o2 




0.51 





f«60000 





150 


YBa2Cu306.85 


Er,Yb EPR0 


9460 


0.2 




25 


160 
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